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We study the application of the coherent-state path integral as a numerical tool for 
wave-packet propagation. The numerical evaluation of path integrals is reduced to 
a matrix-vector multiplication scheme. Together with a split-operator technique 
we apply our method to a time-dependent double-well potential. 



1 Introduction and Definitions 

Throughout this paper, we will consider a standard Hamiltonian 

H = f + V=^ + V(Q) (1) 

for a system with one degree of freedom, described by the momentum operator 
P and the position operator Q. 

A coherent state \a) may be defined by means of harmonic oscillator 
creation and annihilation operators a) and d, respectively, 

' (^/^Q + i—^P), (2) 



^/2h " ' V TOW o 

through 

| a) := exp(aa^ — a*o)|0) , a € C, (3) 

where |0) is the normalised ground state of the harmonic oscillator Ho = 
hojQ [cva + 1/2) , and the exponential is a displacement operator. (Note that 
the frequency u)q, and hence the characteristic length scale (h/mmo)^ 1 ^ 2 , is 
completely arbitrary and can be used as an adjustable parameter.) 
The time evolution of a coherent state under Hq is simple: 

e -^ A °\a) =c- l -^\ae- [uJot } . (4) 

For an operator A we define the antinormal symbol A— (a) implicitly by 
the relation (d 2 a := d Re a dim a) 

A= [ f^\a){a\A.(a). (5) 
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By virtue of the generalised Trotter formul 

lim (F(t/N)) =c-' ltH ' h , (6) 

where F(t) is any operator-valued function with the two properties 
F(t = 0) = 1 

F(t = 0) := lim j (P(t) - l) = -iH/h (7) 

we are able to define the antinormal coherent-state path integral (ACSPI). 
However, to exploit the trivial time development of a coherent state under a 
harmonic oscillator (eq. |]), we consider F of a generalised split-operator type: 

F(t) := e"^ G{t) e~^^° (8) 

with an operator G such that (f7j) holds. We represent G by eq. (||) and define 
the ACSPI 

(a\e- itfl/n \a') :=e' iuJ0t/2 - (9) 
r m j2 N N 

iim / ^...^r[K e+i " 0t/2 l^+i e " iW0t/2 )n G -K;^) 

N^oo I 7T 7T - LJ - - LJ - 

a = a e +i " ot / 2 ,a' = ajv+ie -i^t/a. 
Here we use 

K l—\t\ n 

G(t)--=T,-^r^ 1 /nr (io) 

n=0 

2 ~ 

where Hi = V(Q) — ^rp-Q 2 is the anharmonic part of the Hamiltonian H and 
Ke {4,5,. ..,10}. * 

An analogous normal coherent-state path integral may also be defined] 2 ] 
but will not be considered here. 



2 Numerical evaluation 

For numerical evaluation of the ACSPI (eq. ^) we have to stop at a finite 
(Trotter-)number N of integrations and we perform each integration by a 
quadrature formula 
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Figure 1. The Q-expectation value as a function of time for a particle in the potential of 
equation (ha) is shown. In absence of an external field (S = 0) the particle is tunnel- 
ing through the barrier due to nearly degenerate energy eigenvalues. The destruction of 
tunneling for S ^ is due to degenerate quasi-energies in the Floquet picture.fj 



with fixed sets of abscissas otj and weights Wj > 0, and defining a matrix 
P(t/N) and a vector v by their elements 

Pij = V^e"^ {a lC +i ^\a 3 ) G_ (a,e +i & ; t/N), (12) 

and 

Vj := y/wjiajli)), (13) 

the discretized version of a single time step propagation becomes a matrix- 
vector-multiplication 

v[ := ^(a^e-^W ~ £ P tj (t/N) Vj . (14) 

3 

The vector v' represents the wave packet propagated by the time t/N. 

For time-dependent potentials we choose the time step small enough to 
treat the potential as time-independent during the time t/N. 
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3 Application to tunneling phenomena 

We apply our method to a symmetric double well potential with an external 
time-periodic linear potential: 



The time-independent case (S = 0) shows the phenomena of tunneling (sec 
fig. 0): A wave packet starting in one welLmoves through the barrier into 
the other well. However, Grofimann et al.B observed that application of a 
time-periodic linear potential with the right strength S and frequency to can 
localize the particle in one well. 

In figure [j] we show the ^-expectation value of the wave packet as a func- 
tion of time. Without external field the particle needs about 2500 elementary 
(single well) oscillation periods to tunnel from one well to the other. How- 
ever, application of a field with strength S = 0.0031 and frequency to = 0.01 
suppresses the tunneling process. It is remarkable that this problem has three 
different time scales: The vibrational period T = 2tt/u>q around one mini- 
mum, the period T = 2007T /u>o of the driving force, and the tunneling period 
T w 2- 10 5 /o;o. This shows the capability of our method to give reliable results 
for problems with time scales ranging over several orders of magnitude. 
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